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ABSTRACT 

Three completely different lines of work have recently led to the conclusion that the 
magnetorotational instability (MRI) may create a hoop-stress that collimates jets. 
One argument is based upon consideration that magnetohydrodynamic turbulence, in 
general, and turbulence driven by the MRI, in particular, is more nearly viscoeiastic 
than it is viscous. Another argument is based upon the dispersion relation for the MRI 
in the context of ID simulations of core collapse. Yet a third argument rests in the 
results of direct numerical magnetohydrodynamic simulations of collapsars and thick 
accretion flows. I elaborate on my previous work regarding the first argument above 
and I briefly discuss how these three sets of results are all related. I also discuss the 
different roles played by the magnetic tension and the magnetic pressure within the 
context of this work. I point out that this leads to consideration of the normal stress 
difference between the hoop stress and the radial stress, in preference to a focus on just 
the hoop stress itself. Additionally, I argue that simulations of thick accretion flows 
and collapsars are not self-consistent if they include a phcnomenological model for an 
MRI-induced viscous stress but disregard these other MRI-induced stress components. 
I comment briefly on the RHESSI observation of polarization in the gamma-ray burst 
GRB0212206. 1 argue that this polarization is consistent with a tangled field, and does 
not require a large-scale organized field. Finally I suggest that the role of magnetic 
fields in creating jets, as described here, should be understood not to work within the 
confines of magnetocentrifugal models of jets, but rather as an alternative to them. 
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1 INTRODUCTION 

It has recently been suggested by a number of authors that 
a magnetic field self-consistently generated by the magne- 
torotational instability (MRI) can collimate jets. This has 
been pointed out, independently, following three completely 
separate lines of work. Two are based in simple analytical 
arguments, and a third rests in the interpretation of numer- 
ical simulation. 

These arguments have been placed in a variety of con- 
texts, including supernovae and collapsars, active galactic 
nuclei, and protostellar jets. Certain aspects of this previ- 
ous work are admittedly contextually specific. None the less 
I argue that MRI-generated hoop-stresses may be a nearly 
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universal aspect of astrophysical jet creation and collimation 
in accretion or collapse scenarios. I explore this possibility 
below, where I highlight similarities in these three lines of 
reasoning. However, I focus primarily below on expounding 
on my previous work, as this has not appeared to date in 
peer-reviewed form. 

In Williams (2001) (hereafter W01) I pointed out that, 
in contrast with purely hydrodynamic turbulence, mag- 
netohydrodynamic (MHD) turbulence in an azimuthally 
shearing environment creates a hoop-stress, 1 in analogy to 
the hoop-stress in viscoeiastic media undergoing azimuthal 
shear. I estimated the magnitude of this stress in a thick ac- 
cretion disk (or flow) using simple dimensional arguments, 
and I suggested that this hoop-stress could collimate jets. 
Although I discussed stresses in the general terms of MHD 

1 More precisely, I pointed out that it creates a positive difference 
between the azimuthal hoop stress and the radial stress. 
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turbulence drawing its energy from azimuthal shear, I specif- 
ically pointed out in W01 (viz, p. 2, paragraph 1) that, as 
should be expected in an accretion scenario, this includes 
the Balbus-Hawley instability, i.e. the MRI. (It is not clear 
what other MHD instabilities exist that may draw their en- 
ergy from the differential shear in the context of accretion 
flows.) This specific point regarding the role of the MRI was 
later emphasized in Williams (2003) (hereafter W03). The 
context of W01 was a generic thick, steady accretion disk 
or flow, and in W03 I discussed protostellar jets, although 
most of the latter argument is broadly applicable to other 
jet phenomena. 

It was suggested by Akiyama et al. (2003), see also 
Akiyama & Wheeler (2003), that a toroidal field generated 
by the MRI could contribute to jet creation in core-collapse 
supernovae. This argument was based upon the dispersion 
relation for the MRI. The MRI is a weak-field instability 
in the sense that the MRI is unstable so long as the mag- 
netic field is less than some critical field strength. In addi- 
tion, the growth rate of the MRI is comparable to the local 
angular velocity fi, which enables several e- foldings of the 
field strength in the course of core-collapse. It was thus hy- 
pothesized that the MRI would amplify a small seed field 
until it grew to equal the critical field strength, and it was 
argued that the resultant magnetic field would contribute 
to jet creation. This tentative conclusion was bolstered by 
ID core-collapse simulations in which the MRI was treated 
heuristically through its dispersion relation. 

Collimation by an MRI-generated field has also been ar- 
gued to be happening in several recent direct numerical sim- 
ulations (DNS) of both transient and quasi-steady accretion 
flows. This is particularly intriguing because DNS simula- 
tions, by definition, make no assumptions about the MRI 
or the ensuing turbulence. The first simulations in 3D seem 
to show collimation in the form of a funnel surrounded by 
high magnetic pressure regions (Hawley et al. 2001, Hawley 
& Balbus 2002). This is quite different from the collimation 
by hoop-stress that I have argued. Meridional plane (2.5 D) 
simulations of accretion flows in a variety of contexts seem 
to show collimation by stresses that the authors argue are 
generated by the MRI (Kudoh et al. 2002, Proga & Begel- 
man 2003, Proga et al. 2003). Further numerical simulation 
and study should help clarify the role of the MRI in this 
collimation. 



2 COLLIMATION BY TURBULENCE 
2.1 A preliminary zeroth route 

Tables for the components of the Reynolds stress and turbu- 
lent Maxwell stress tensors in shearing-sheet simulations of 
the MRI have been provided by Brandenburg et al. (1995), 
Hawley et al. (1995) and Hawley et al. (1996) (hereafter 
BNST95, HGB95 and HGB96). The stresses are due largely 
to a tangled magnetic field, and the streamwise stress, cor- 
responding to a toroidal hoop stress, is positive (indicating 
tension) and larger than the (off-diagonal) viscous stress; 
in contrast, the radial normal stress is negative (indicating 
pressure). Tables for the force due to a stress tensor (i.e. 
its divergence) can be found in a variety of standard ref- 
erence sources: A positive hoop-stress Wee and a negative 



radial stress Wrr result in a cylindrically radially-inwards 
(collimating) force density F R 2 ' = —(Wee — Wrr)/R. This 
leads immediately to a zeroth conclusion that the MRI can 
create a collimating force 2 , and that it can do so through a 
turbulent, tangled field, by direct inspection of the results 
of BNST95 and HGB96. The primary serious potential ob- 
jection to this is that the magnitude of the collimating force 
may not be significant, since the MRI saturates at a rather 
high plasma /3 P in these simulations. There remain several 
additional questions, such as how magnetic pressure affects 
this collimation. In particular note that the radial stress 
Wrr is due largely to magnetic pressure, and the gradient 
term OrWrr in the divergence of the stress [F^\ infra] 
in a real accretion flow can not be determined from local 
shearing-sheet simulations. The remainder of this section 
discusses these and other observations and questions regard- 
ing the MRI in greater detail, with the help of some simple 
turbulence modeling considerations, following and extending 
the discussion in W01 and W03. 



2.2 Some simple turbulence modeling 
considerations 

2.2.1 Introduction and notation 

By viscosity, unless otherwise noted, I mean shear viscosity 
and not bulk viscosity. For the purposes of this paper, by 
viscous stress I mean that part of a stress tensor w, turbulent 
or otherwise, which may be written in the explicit tensor 
form 



X7v+(X7v) T - ^(V-v)/ 



(1) 



where v is the velocity and / is the identity. The only cases 
considered here are rectilinear shear, such as v oc yx in 
Cartesian coordinates, and azimuthal shear such as in an ac- 
cretion disk with Rfl 3> vr. In these two cases, the viscous 
stress tensor above reduces to a single, off-diagonal compo- 
nent in the respective coordinate systems, namely w xy and 
WRe- It is assumed that this stress component is due entirely 
to viscosity, so that the value of the kinematic viscosity be 
it a molecular or a turbulent viscosity, is given by v — w/7 
where 7 is the shear rate and w is the appropriate stress 
component above, and care is taken so that the sign of v is 
consistent with the tensor equation (1) above. 

My approach here as elsewhere is to ignore the mean 
field entirely, as counterpoint to analyses that ignore the 
turbulent field entirely. Thus unless noted otherwise I as- 
sume throughout the remainder of Section 2 that the mean 
field (B) is identically zero, and the magnetic field exerts a 
dynamical influence entirely through a tangled field B' . I de- 
note the averaged Faraday tension as AnMij = (B^B'j), and 



the turbulent Maxwell stress M, is given by Mi 



Mi. 



(l/2)MkkSij- The Reynolds stress is Rij = (pv'iVj), and the 
full turbulent stress tensor W is given by Wij = —Rij+Mij. 
Note that it is commonplace to refer to the funicular 



2 When discussing collimation here and in my previous work, it 
should be clear that I am not discussing far-downstream colli- 
mation, such as in traditional magnetocentrifugal mechanisms. 
Rather, I mean collimation in the jet-formation region, including 
the equatorial plane of accretion. 
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turbulent stress My as the Maxwell stress. This is strictly 
speaking incorrect; one must not neglect the contribution of 
the magnetic pressure term to the turbulent stress tensor. 
Inclusion of the turbulent magnetic pressure term does not 
affect the viscous term nor the normal stress difference, but 
it is none the less important in a full consideration of the 
dynamics of how MHD turbulence collimates jets. 

Also for future reference, distinguish here between the 
angular velocity or Coriolis rate fi, the shear rate 2 A = 
—ROrQ (here A is the first galactic Oort constant), and the 
vorticity (1/ R)d R (R 2 Q.) = 2(Q - A). In the case fi oc R~ q , 
these may be written fl, qQ, and (2 — q)Q, respectively. 

2.2.2 Failure of a -prescription and existence of normal 
stress difference 

The original way (Shakura & Sunyaev 1973) to write the 
Shakura-Sunyaev viscosity prescription, under the assump- 
tion that vr <C v$, is that the viscous component of the 
turbulent stress tensor is proportional to the density times 
the sound speed squared, Wrq = —a ss pc 2 . This represents 
only one of the six independent components of the turbulent 
stress. Through vertical integration, for a quasi-2D (R-6) 
disk, these six reduce to three, namely the aforementioned 
viscous stress and the on-diagonal stresses Wrr and Wee- 
These on-diagonal components of the turbulent stress are 
usually ignored. If anything, it is typically assumed that 
they consist either simply of the so-called turbulent pres- 
sure, Wrr — Wee = — Pturb, or are due to a turbulent bulk 
viscosity. Either assumption predicts that the on-diagonal 
stresses are equal to one another. Particularly in the case of 
MHD turbulence, this may be highly inaccurate. 

On-diagonal stress may be dynamically significant even 
when the turbulent pressure and bulk viscosity are not, be- 
cause of normal stress differences. It is therefore useful to 
distinguish between what I call the primitive a-viscosity 
prescription above where only the viscous stress Wrb is mod- 
eled, and the extended a-viscosity prescription in which the 
full turbulent stress tensor is modeled as a purely viscous 
stress, plus perhaps a pressure term. 

The shearing-sheet simulations of the MRI mentioned 
above clearly show that: (i) the stress is dominated by the 
turbulent Maxwell stress, rather than the Reynolds stress, 
and (ii) as pointed out above, the cross-stream stress and 
the streamwise stress are of opposite sign, and both are ac- 
tually larger in magnitude than the viscous stress, creating 
a significant normal stress difference, in gross contradiction 
of the extended a-viscosity prescription. These two facts are 
connected (Williams 2004b). 

It should be clear by examining BNST95, HGB95, 
HGB96, as well as Matsumoto & Tajima (1995), that the 
magnetic field that is being produced by the MRI is being 
dragged by bulk shear of the fluid, as ideal MHD tells us it 
should, and it is statistically aligning itself with the direction 
of mean shear. It is precisely by such dragging of field lines 
that the Maxwell stress in a turbulent medium produces a 
quasi-viscous off-diagonal stress, but why should the drag- 
ging of field lines stop there? Let us imagine, as I did in W01, 
a process where a turbulent or tangled field is constantly be- 
ing created, distorted by shear, and ultimately dissipated, as 
appears to be the case in these shearing- sheet simulations 
(see fig. 2 in Williams (2004b)). Note that in equilibrium, 



when the turbulence is saturated, this sequence is a logical 
sequence, not a chronological sequence. 

In fact, compare the magnetic stress My from eq. (4) 
of W01 to the normalized stress My of eq. (24) of BNST95, 
keeping in mind the switch x y in going between the two 
papers, as well as the change of sign of M X y according to 
the differing orientation of the shear. A least-squares mini- 
mization of the difference between each of the four nonzero 
members of the six independent stress components taken in 
turn gives a fit of the normalized My of W01 to the My 
results of BNST95, written in the coordinate system of the 
latter, which corresponds to pure azimuthal shear in a disk 
with ordering (R,0, z), of 



(My )woi — 



(My 



\ norm 
ij JBNST95 




(2a) 



(2b) 



for a value of the normed relaxation time 7s = 6.9 in the 
notation of W01. 

In the case of the MRI the sequence of creation, dis- 
tortion, and dissipation mentioned above is part of a larger 
feedback loop where the shear-aligned magnetic field is un- 
stable to the MRI (locally, through the combined action of 
the shear operator and the Coriolis operator) , creating more 
turbulence and completing the feedback loop: 

(3) 



(Mee+)MnR 




shear -\-rotn (MRI) 



Mee{-M RR ) 



In principle, the feedback loop may be closed by other 
sources of turbulence, such as convective or inertial insta- 
bilities, or by MHD instabilities other than the MRI; it is 
in the nature of turbulence to stretch and contort material 
lines. However, note that passive turbulent dynamo action 
(such as magnetoconvection) appears to be much less pow- 
erful than the active turbulent dynamo action of the MRI 
(see HGB96 regarding this point). 

Indeed, let it be clear that the process I have described 
here and previously (W01, W03, Williams (2004a) Williams 
(2004b)) depends upon dynamo action in the sense of a tur- 
bulent dynamo as described by HGB96. Misunderstanding 
on my part regarding the definition of the word dynamo, 
a word which appears often to be taken to be synonymous 
with processes that can create self-sustaining large-scale or- 
dered fields, led me to state otherwise in Williams (2004a), 
which was incorrect. On the other hand, the process does 
not depend upon a mean-field dynamo. 

In fact, a large-scale magnetic field will complicate 
this picture, as will other symmetry-breaking terms such 
as VP x f2. None of these are present in HGB96, nor are 
they present in the analysis here. A mean magnetic field is 
present in HGB95; the orientation of this field seems to af- 
fect the effective relaxation time for the turbulence, as I note 
in W03; more significantly, the presence of this field affects 
the saturation of the MRI. Ignoring the mean field and other 
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symmetry-breaking terms for the present purposes however, 
I find that simple viscoelastic models of the stress produce 
a stress tensor W that I write here in terms of the ther- 
modynamic pressure P for a disk in which v R <C vg as (cf. 
Williams (2005) infra eq. 8) 

W/P = -a I -a x R0 + a 2 [00 - (1/3)7] (4) 

where the set {ai} consists of dimensionless constants, and 
R and are unit vectors. That is, the stress is the sum of 
a turbulent pressure equal to and parameterized by ao, a 
turbulent viscosity parameterized by on ~ o ss (given that 
P ~ PCs), and a streamwise normal stress difference, 3 4 
due to the turbulent elasticity, parameterized by 02. The 
form chosen for this final term simply makes it traceless in 
3D. For comparison with Williams (2005), ao = Pturb/P, 
Qi = Mturb(- Rd R il)/P, and o 2 = 2( turh (Rd R n) 2 /P. For 
future reference observe that (i^ 1 = (3/2)ao(— Mkk/Wkk) = 
3ao{A4kk/Wkk) — (1.5 to 2.3) ao ~ ao. The lower number 
(1.5) is derived from HGB96 and the higher number (2.3) 
from BNST95. 

Note that the model for the MRI turbulent stresses 
given by Ogilvie (2003) produces a different form for the 
equilibrium turbulent stress, because of the explicit inclu- 
sion of the Coriolis effect through a Coriolis operator in his 
eq. (23). This Coriolis operator rotates the Reynolds stress 
about the vector tt in the local corotating reference frame, 
and it is most likely a crucial ingredient in a complete local 
dynamical model for MRI driven-turbulence. However, this 
difference does not materially affect the conclusions here: 
The existence of significant normal stress differences is an 
unavoidable feature of any realistic turbulence model for the 
MRI in most, if not all, contexts, and this gross feature is 
present in the model of Ogilvie (2003) as well. 

Let me be the first to point out that there are three 
variables to model (the three stress components) , and I have 
done so using three parameters, so eq. (4) not predict any- 
thing by itself, although note this prescription works quite 
well for describing the full six-component stress tensor in 3D 
as well ( W03) . The reason this is a meaningful expression is 
that a2 > ai, and that it makes physical sense to separate 
the normal stress into an isotropic and an anisotropic part. 

2.2.3 Magnitude of the normal stress difference and the 
effective Weissenberg-Deborah number 

The ratio of the normal stress difference to the viscous stress 
(02/01 in eq. 4) is roughly proportional to the ratio of the 
effective stress relaxation time s to the shear time scale 
(2 A) -1 , with the exact constant of proportionality being a 
model-dependent quantity. For example, the Maxwell model 
of Ogilvie (2001), which I write here as 

s{i> t W) + W = u[Vv + X7v T ] + (vb - ^uj (V ■ v)I (5) 

3 The first normal stress difference in the case of pure azimuthal 
shear is N\ = Wgg — W RR . The second normal stress difference is 
N2 = Wgg — W zz . I am only concerned with JVi here; henceforth 
the first normal stress difference will simply be called the normal 
stress difference. 

4 Normal stress differences are not coordinate-invariant, but they 
are none the less useful and physical quantities within the context 
of the work here. 



predicts {Wgg - W RR )/W R g = 2s(2A). (Note that V t is the 
tensor generalization of the vector advection operator for 
the magnetic field B, see Ogilvie (2001), W01, et seq.) As 
well, the simple a-S model (model B) of W03, namely 

s(p t M) + M = aI (6) 

predicts (Wgg — W RR )/W R g = s(2A), under the approxima- 
tion that 

Wee-W RR ~Mee-M RR = Meg-MRR. (7) 

The a-5 model was not explicitly given in W01, although the 
construction in W01 gives a stress of exactly the same form, 
in the case of steady shear that is appropriate here. The ap- 
proximation in eq. (7) follows from the assumption that the 
turbulent Maxwell stress normal difference is much larger 
than the Reynolds stress normal difference, as a result of ei- 
ther the quicker relaxation to isotropy of the Reynolds stress 
compared with the turbulent Maxwell stress, or (more signif- 
icantly) the assumption that the turbulent magnetic energy 
is greater than the turbulent kinetic energy. This second 
assumption is important to the model of Ogilvie (2001) as 
well; otherwise the use of the magnetic tensor advection op- 
erator for the advection of the full turbulent stress would be 
unjustified. 

Independent of the question of the value of the viscosity 
in a disk, then, is the question of the value of the ratio of 
the normal stress difference to the viscous stress. Henceforth 
in this paper, I define the effective Weissenberg-Deborah 5 
number We to be equal to the ratio (Wgg — Wrr)/\Wr$\. 
In the notation of eq. (4) above and of Williams (2004b), 
for reference, We = 02/01 = 2(£//i)(2A). It is also useful to 
define no = 00 / ai . The relative magnitude of the members 
of the set {o^} are then determined by 7ro and We. 

Basic dimensional arguments suggest that We is of the 
order of unity: take the limit of an inviscid, perfectly con- 
ducting fluid, with zero mean field. Under the assumption 
that the Coriolis parameter f2 and shear rate are compara- 
ble, the only local background dimensionful quantities avail- 
able to construct a time scale are f2 _1 and, given a length 
scale L (such as a scale height), a sound-crossing time L/c s . 
For a thin disk these two are also both comparable, and 
this suggest that, barring some large or small dimensionless 
factor (such as the magnetic Prandtl number), the effec- 
tive relaxation time is comparable to the shear time scale 
(Ogilvie 2001, W01, Ogilvie 2003, W03, Williams 2004b). 
One should then expect the streamwise normal stress differ- 
ence to be comparable to the viscous stress, in gross con- 
trast with the extended o- viscosity prescription, and this is 
what is seen. On the other hand, note that almost the exact 
same line of reasoning has been used, historically, to argue 
that the original Shakura-Sunyaev o ss ~ 1, whereas the ac- 
tual value of this parameteris typically orders of magnitude 
smaller. For reference, note that the two time scales above 
start out equal in the simulations of HGB95; at the end of 
the simulations, c s /L ~ 5f2 for the vertical field fiducial run 
and c s /L ~ 4fi/3 for the azimuthal field fiducial run. 

5 Note that rheologists often distinguish between the Weis- 
senberg number and the Deborah number; in principle it is pos- 
sible to have one large and the other small, see Section 4.2 of 
Phan-Thien (2002). Here I make no such distinction. 
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A much stronger argument can be made by appealing to 
the action of the MRI as a feedback loop, as in diagram (3). 
This feedback loop can only operate if the relaxation time 
of the process is long enough to create a substantial dif- 
ference between the streamwise stress Mgg and the cross- 
stream stress M RR ; otherwise the feedback loop will not 
close. This suggests that We may be arbitrarily larger than 
unity but it is hard to see how it could be smaller than unity. 

Even more significantly, an analysis of BNST95, 
HGB95, and HGB96 shows that that 02 is not just compara- 
ble to but in fact larger than both ao and ai . Some analysis 
of HGB95 is found in W03. For HGB96, the full stress tensor 
(taking into account the often neglected contribution of the 
magnetic pressure to the turbulent Maxwell stress, supra) 
from a naive average of their Rl, R2, R3, R4, R6, and R7, is 
(W RR ,W B B,W zz ,W m ) = (-0.038,0.012,-0.028,-0.015). 
The remaining two stress components, Wg z and W zR , are 
not provided, although symmetry suggests they are zero. 
I find ao, on, and 02 to be 0.018, 0.015, and 0.046, respec- 
tively; this gives We = 2.98 and tt = 1.19. The model eq. (4) 
is also an excellent fit to the full stress tensor provided in 
BNST95 as well. The full normalized stress tensor (taken 
from their eq. (24), which provides all six components of the 
separate Reynolds stress and the Faraday tension contribu- 
tion to the turbulent Maxwell stress My), again, correcting 
for the neglected contribution of the magnetic pressure to 
the turbulent Maxwell stress, and assuming a ratio of mag- 
netic energy to kinetic energy of 6.67 taken from averaging 
over all runs given in their table 1, and normalized so that 
Tr(W) = -1, is 



/-.780 -.159 \ 
(Wijjw,™ = -.159 .487 (8) 
\ -.707/ 



under the coordinate ordering (R, 9, z) for the columns and 
rows. This gives a value We = 02/01 = 7.7, larger than the 
value of We obtained from HGB96. On the other hand, it 
is possible that the difference in We is simply reflective of 
the differing aspect ratio of these simulations. As well, of 
course, future higher fidelity simulations may produce even 
different results. 

Significantly, the variation in the value of We obtained 
in HGB96 compared with BNST95, i.e. roughly 8 versus 
roughly 3, is much less than the variation between these two 
studies of the quoted value of a ss , which differ by a factor of 
almost 10. This observation is in line with the observations 
of HGB95 [see their discussion in paragraph 1 of p. 749, 
and their eq. (16) and eq. (20)] that, according to them, the 
traditional a parameterization 'does not seem as appropri- 
ate here' because the viscous stress correlates more strongly 
with the magnetic pressure than with the gas pressure. Since 
the field, and hence the magnetic pressure, is dominated by 
the streamwise toroidal field, their observations combined 
with the observations in the paragraph above suggest that 
the ratio of the stress components, and in particular, I would 
argue, the ratios ao : ai : a.2, are a more robust result of 
shearing- sheet simulations than the absolute value of the 
viscous stress as parameterized by the traditional a ss . 



2.3 The form of the stress in a thick accretion 
flow 

2.3.1 Steps towards a model 

In this section and the subsequent two sections, I describe 
basic steps towards a turbulence-based model of jets. A more 
complete description of such a model will be given in a later 
paper. The stresses and forces to be considered here are the 
turbulent stress and the forces due to gravity, the pressure 
gradient, and the centrifugal force. The centrifugal force is 
the result of the 99 component of the inertial stress tensor; 
the remaining terms in this stress (RR, Rz and zz) that ap- 
pear in the meridional equation of motion create a ram pres- 
sure which may create forces in the meridional plane that 
are of comparable magnitude to the other forces described 
here and which may also aid in collimation. However, for the 
sake of simplicity here, these forces are ignored. 

2.3.2 The fundamental assumption 

From the values given above for the set of numbers {cti}, in 
a Keplerian thin disk the resultant collimating force is even 
less than the radial thermal pressure gradient, and thus neg- 
ligible. Were this not the case, the thin-disk approximation 
would fail to be self-consistent. 

Of much more interest is the case of a geometrically 
thick disk or flow. It is not yet clear how turbulence saturates 
in this context. A global accretion flow is different from a lo- 
cal shearing-sheet simulation in many respects. Among other 
things, the restrictive assumptions (quasi-periodic shearing 
box boundary conditions with no handedness) necessary to 
perform the shearing-sheet simulations of the MRI may af- 
fect the plasma j3 p at which the MRI saturates. In fact, as I 
discuss below, preliminary global simulations show the MRI 
to be more powerful (by which I mean j3 p is smaller) than 
the local shearing-sheet simulations do. 

My approach here, as previously, is to fix the value of 
We (and 7ro) for the reasons described above, but to let 
the overall magnitude of the viscous stress be given by a 
free parameter. Ignorance of the behavior of turbulence in a 
thick disk (or flow) is thus subsumed entirely into the vis- 
cosity prescription, and it is assumed that the other stresses 
(streamwise normal stress difference and turbulent pressure) 
occur in proportion to the viscous stress. Previously (W01) 
I fixed We to be of order unity. If instead we take We ~ few 
as above, then so much the better. 

It remains to be seen what is the absolute value of the 
viscous stress W R e and how this compares with other dimen- 
sionally similar quantities. In the next few sections I provide 
an unfortunate but necessary diversion into the various pre- 
scriptions for the viscosity. Largely, the variations in ways of 
writing the viscosity are irrelevant to the gross conclusions. 

2.3.3 The form, of the viscosity in W01 and W03 

The specific scenario discussed in W01 and W03 was a flow 
driven by the viscous spin-down torque on a central object 
completely embedded in a thick accretion flow (I did not use 
the word flow but rather disk, but let it be clear that this 
in no way should be taken to imply equilibrium as in the 
classic studies of thick equilibrium tori). In W01 I chose to 
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write v in the form of a proportionality constant (a) times 
a squared length scale times a rate. The rate I chose was the 
shear rate at the fiducial radius of the central accretor R*, 
so that v oc (2A) t = gf2*. The length scale was somewhat 
idiosyncratically written in W01 as y/]3R t , and j3 was a free 
parameter (not the plasma beta) allowed to vary to take 
into account the uncertainty as to the correct length scale, 
and thus viscosity. (To avoid confusion with the plasma /3 P , 
henceforth label j3 from W01 as /3 W01 ; also label a from 
W01 as a wol . Since a wol is also a free parameter and the 
two only appear in the form of the product c* W01 /3 W01 , this 
second parameter /3 W01 is redundant in the notation of W01. 
Regrettably, this redundancy obfuscated the primary thrust 
of the paper.) In particular, I allowed the length scale (here 
let me call it L) to vary between J?* as a lower bound and the 
radius of the thick disk as an upper bound, which I denoted 
r m R* (note the omitted minus sign in W01 and W03): 

R R 2 ( ^ \ 

awoiP ™ 4 \dbiRj m 

= Q L 2 (^-) W 
W01 \d\nRj t 

= u W01 qL 2 Q.* 

The length scale /3 W01 R, , or L, was assumed for simplicity 
to be a globally defined quantity that does not vary with 
position, and so the viscosity above should be understood 
to be constant throughout the flow, for the sake of argument. 

The reason for the use of the parameter /3 W01 was simply 
that it is sometimes taken to be implicit, when writing the 
viscosity in the form v ~ £ 2 uj, that the length scale t and 
the rate lo are in some dynamical sense characteristic of the 
turbulence. I wished to emphasize that, for a central object 
completely embedded in an accretion disk that is potentially 
much larger than the central object, this length scale could 
be much larger than the radius of the central accretor, but 
no larger than some effective characteristic length scale of 
the thick accretion region. 

This suggests, if a wol in eq. (9) is of the order unity, 
that the product ct W01 f3 W01 in principle may be much larger 
than unity. Note that this does not necessarily imply any- 
thing about the value of a ss , which may still be of the order 
of or much less than unity, see below. 

2.3.4 Shear versus vorticity versus Coriolis parameter 

The form I adopted above posits the viscosity being pro- 
portional to a shear. Abramowicz et al. (1996) performed 
shearing-sheet simulations of the MRI with variable effective 
q and found that the viscosity is proportional to the ratio of 
shear to vorticity, Q./2A, or equivalently, v oc q/(2 — q)c 3 H. 
This is an important consideration to take into account for 
a model of the MRI for future work on thick disks and flows. 
However in the rough work presented here this only intro- 
duces a factor of order unity which is negligible relative to 
the other uncertainties and approximations in this work. 

2.3.5 Global versus local viscosity prescriptions 

Historically the turbulent viscosity has been variously pa- 
rameterized in terms of purely local quantities, purely global 
quantities, and mixtures of the two. For example, the 



Shakura-Sunyaev viscosity parameterization in its original 
form, supra, is a local definition, whereas the viscosity pa- 
rameterization of Lynden-Bell & Pringle (1974), namely v oc 
£l,R 2 , is in global form. Although Lynden-Bell & Pringle 
(1974) do not define an alpha, it is convenient to do so us- 
ing their formulation, so let me here define v = a^pfl^R 2 . 
From this one may define an effective Reynolds number 
R-e = 1/cw. 

Popular mixed forms for writing the viscosity prescrip- 
tion include v — aH 2 Q, and v — ac s H, where H is some 
measure of the disk thickness. These definitions are more or 
less equivalent to the original Shakura-Sunyaev prescription, 
in the case of a thin disk. This redundancy disappears in the 
case of a thick disk or flow, however: it is not clear what is 
the natural generalization of H and Q in a thick accretion 
flow, because there are are several nearly equal length scales 
and rates that are habitually conflated in thin-disk theory 
into these two quantities. For example, H in thin disk theory 
is neither strictly local nor is it global, it lives somewhere in 
between. It is variously defined to be a representative pres- 
sure scale height, a density scale height, the second moment 
of the vertical density distribution, or the half-thickness of 
the optical depth t — 2/3 surface. It is not sufficient to gen- 
eralize it to a local quantity by defining it to be a pressure 
scale height, say, because strictly speaking the local pres- 
sure scale height in a thin disk is infinite at the midplane. 
Similarly, Q, in thin disk theory is vaguely synonymous with 
the orbital frequency, the Coriolis parameter, the Keplerian 
orbital frequency, the shear rate, and the vorticity; none of 
these quantities are necessarily nearly equal one another in 
a thick disk. This opens up the possibility that a may be 
anomalously larger or smaller than the accepted value, de- 
pending upon what one means by a, and what part of a 
three-dimensional accretion flow one is examining. 

The choice in W01 was to write the viscosity in a purely 
global form, similar to the Lynden-Bell-Pringle formulation 
above. In particular, a^p = qct W01 (3 W01 . Henceforth, I dis- 
pense with Q wol and /3 W01 in preference for a mp . 

The values of q ss and a mP are connected in a non- 
trivial way. Under the approximation that P ~ pc 2 s and 
2 A ~ Q., one finds 

«ss " («?). /n. 1 ' 

where vg = RQ. Near the surface of the central accretor, 
then, Qujp/Qss — c 2 /v 2 . For R> R*, one may have a^p 3> 
a ss if c'i > v%. 

Having now discussed the variations in writing the vis- 
cosity, I now return to the physics of collimation. Follow- 
ing W01, for concreteness, assume that in the thick disk, 
fl ~ R~ q , and in particular f2 may differ greatly from the 
Keplerian value Qk ■ 

2.4 The ordering of forces and speeds 

The force density that arise from taking the divergence of 
the stress tensor may be written as the sum of three forces, 
see the three terms on the right hand side of eq. (21) below, 
respectively called F"\ F^ 2 \ and here. In this paper I 
do not consider F^. 

The normal stress differences create a cylindri- 
cally radially inwards specific force f^(~ Fr^/p) = 
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-(We)W Re /R 



-(We)v\d R n\ (c.f. W01, eq. 6). In the force is 



(2) 



Shakura-Sunyaev local form, this may be written f R 
— (We)a ss Cg/R, and in the Lynden-Bell-Pringle global 
form, f R 2) = -(We)a tBP ^n«(n/_R). For We ~ 1 and for 
Re < 1 (i.e. a mP > 1) as assumed in W01, the collimating 
force due to the MHD turbulence is clearly significant. 

Consider here the relative strength of various specific 
forces. For simplicity here, assume that our test fluid element 
is close to the equatorial plane, z/R <C 1. The set of forces 
I wish to consider here consists of the force due to the first 
normal stress difference f R \ the centrifugal force f c = Q 2 R, 
and the force due to gravity f g = ii%R. The force due to the 
turbulent pressure (which is essentially the same as the 
force due to the magnetic pressure because Mkk 3> Rkk), 
and the force resulting from pressure fp = — (8rP)/p will 
be discussed in the subsequent section below. 

Before proceeding with the evaluation of forces, how- 
ever, it is handy to describe the collimation stress in terms 
of a wavespeed, just as Akiyama et al. (2003) make use of 
the Alfven speed to describe the strength of the magnetic 
field. A tangled field produced by a turbulent conducting 
medium supports Alfven-like waves much like a viscoelas- 
tic medium supports elastic waves (Gruzinov & Diamond 
(1996), Schekochihin et al. (2002), Williams (2004b)). The 
speed of these transverse elastic waves n tow is dependent 
upon direction, for an anisotropic stress tensor. The wave 
speed due to the toroidal component of the tangled, turbu- 
lent field is similar to the expression for the Alfven wave 
speed, «tewW = (BgBg) /Anp. Then, for the toroidal mode, 



2 

"tew 



P 



(11) 



Note that 



Mee = ( 2a + -012 



27To + - We ) a ss pc 3 



(12) 



- (We)a ss pc s 
and therefore for the toroidal mode 
vl^O) ~ (We)a s 



(13) 



Typically the shearing-sheet simulations show (We)« ss is 
small, on the order of l/(few), e.g. 1/20 for HGB96. It 
is conceivable that high-fidelity global simulations will pro- 
duce markedly different results; perhaps it is possible that 
(We)a ss becomes larger than unity. Incidentally, note that 
the principle axes of the stress tensor Mtj do not coincide 
with the coordinate axes, so an initially prograde or retro- 
grade azimuthal wave will be refracted radially inwards or 
outwards respectively. Henceforth let it be understood in 
this paper that by »te W I mean the speed of the azimuthally 
traveling wave Utew(S). Also note that i>t ew is equal to the 
energy per unit mass stored in the azimuthal part of the 
turbulent magnetic field. The relative magnitude of the var- 
ious forces can be stated in terms of the relative magnitude 
of various speeds, or the relative specific energy densities as 
well, as given below. 

(2) 

The condition for the turbulent collimating force f R to 
dominate (i.e. have a larger magnitude than) the centrifugal 



%<(We)a ss , 



or, equivalently, 



R 



(14) 



(15) 



Even more succinctly, eq. (14) may be restated with the 
help of the wave speed expression above, so the condition 
becomes 



(16) 



That is, the radial force due to the turbulent normal 
stress difference dominates the centrifugal force when the 
azimuthal flow velocity becomes sub-(quasi)-Alfvenic, or 
equivalently, when the energy in the shear-aligned magnetic 
field is greater than the kinetic energy of rotation. Inter- 
estingly, this critical field strength is the same as the MRI 
saturation field given in eq. (8) of Akiyama et al. (2003). 

(2) 

The condition for the turbulent collimating force f R 
to dominate gravity, in terms of the Keplerian orbital speed 
v K , is 



4 



< (We)a s 



or, equivalently, 



TT < ( We K 



(17) 



(18) 



Again, more succinctly, eq. (17) may be restated with the 
help of the wave speed expression above, so the condition 
becomes 

GM» 



R 



2,2 

V K < "tew 



(19) 



That is, the radial force due to the turbulent normal stress 
difference is stronger than gravity (but pushes in the same 
direction!) when the energy per unit mass stored in the tur- 
bulent azimuthal field exceeds the gravitational binding en- 
ergy per unit mass. 

Of course, there is nothing magical at all in eq. (16) and 
eq. (19). The question is whether the turbulence is ever so 
strong that these conditions are ever approached. Roughly, 
for collimation by turbulence to work, it appears likely that 
we must have that v 2 c „ is a substantial fraction of vg, and 
preferably 

v'g < "tow j 

and then if we assume (We)a ss < 1 this implies 



vg < < <£. 



(20) 



The end result, perhaps not surprisingly, is that the 
turbulent elastic radial force is unequivocally impor- 
tant when the rotation and the Keplerian orbital speed be- 
come sub- (quasi)- Alfvenic as in eq. (16), and furthermore 
subsonic as well if we assume as is usual that c 2 > v 2 ew . 
These are not conditions that are ordinarily met in a thin 
disk, but it is conceivable that they are met in portions of 
a thick disk, where radial pressure support becomes signif- 
icant. Note that the Shakura-Sunyaev formulation implies 
that turbulent pressure and hence magnetic pressure occurs 
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in proportion to the thermal pressure. Even if this is not 
precisely true in reality, it suggests that if we are to dis- 
cuss radial thermal pressure support, then we should as well 
address radial magnetic pressure support. 



2.5 The role of magnetic and turbulent pressure 

This brings me now to address one final potential objec- 
tion to collimation by turbulence-induced magnetic hoop- 
stresses, namely that magnetic pressure might push out and 
counteract the inward pull of the hoop-stress. 

This was explored in W01. As I describe in further de- 
tail here, the existence of magnetic pressure may actually 
be beneficial to the jet-producing scenario I have described. 
It does, however, lead to a focus on the normal stress differ- 
ences, rather than the hoop-stresses themselves. 

The phrases hoop-stress, azimuthal stress and toroidal 
stress are all succinct, equivalent and physically intuitive. 
The phrase 'positive normal stress difference between the 
azimuthal stress and the radial stress' is more precise, but 
not succinct. One may have toroidal hoop-stresses without 
normal stress differences, but this is physically equivalent to 
a negative pressure. Collimation in such a scenario can only 
occur through the radial gradient term below. 

The existence of a positive (tension) turbulent mag- 
netic hoop-stress requires a normal stress difference: the on- 
diagonal turbulent magnetic stress can not be of the form of 
a negative pressure; the normal stresses can not all be pos- 
itive, since their sum is negative semi-definite. Recall that 
the i?-component of the divergence of a symmetrical (stress) 
tensor in cylindrical coordinates is given by 

(v ■ w) R = OrWrr - Wee ~ Wrr + d z w Rz . (21) 

For a purely azimuthal field, spatially stochastic or other- 
wise, the stress Wrr is negative and due entirely to magnetic 
pressure, and it appears in the force equation through its 
gradient. In contrast, the normal stress difference Wee — Wrr 
is due entirely to the magnetic Faraday term, i.e. the ten- 
sion of field lines, and it appears in the force equation only 
through the geometrical factor 1/R, not through its gra- 
dient. In the more general case of a tangled field which 
is predominantly but not entirely azimuthal, not only will 
Wee have some contributions from the magnetic pressure, 
but Wrr will have some positive contributions from the 
tension term as well. However, the normal stress difference 
Wee — Wrr, absent Reynolds stresses, is still due entirely to 
the magnetic tension term. 

Now let us examine outflow collimation in engineering. 
Consider jets from nozzles, such as a garden hose nozzle or 
a rocket nozzle. These jets are confined and collimated at 
their base by the anisotropic elastic stress in the collimat- 
ing medium (i.e. the nozzle). Similarly, the explosive gases 
that propel the slug of a rifle or cannon are effective at pro- 
pelling the projectile becuase of the anisotropic elastic stress 
in the surrounding barrel. The metal in the barrel is radially 
in compression (negative stress) and azimuthally in tension 
(positive stress) . In the case of both the steady dynamics of 
nozzles and the transient dynamics of guns, the collimation 
is thus due to a normal stress difference in the surrounding 
material. In all cases, the anisotropic stress of the surround- 



ing medium reshapes the effect of an isotropic stress (namely 
the pressure) to create an anisotropic outflow. 

On the other hand, it is well to note that underwa- 
ter detonations of explosives may produce jet-like columns 
of water; as well, large thermonuclear explosions in the 
stratified atmosphere of the Earth can cause a rapid ver- 
tical acceleration of material. An analogous phenomenon 
is sup ernovae- fueled outflow in stratified galactic atmo- 
spheres, known as the galactic fountain. Such processes cre- 
ate anisotropic outflows due to collimation by pure pressure. 

The magnitude of the force density due to the gradient 
of turbulent pressure (including magnetic pressure) F^' = 
9r(Wrr) should be of the order of Wrr/R. It remains to 
be seen whether this force is collimating or anti-collimating. 

I assumed in W01 that pressure increases as we near 
the central object, so that the pressure is anti-collimating. 
Then, F (1) and F {2) act oppositely. For the net force to 
be collimating, one requires \F 2 \ > |.F^|. I argued that 
work done by the magnetic curvature forces - i. e. work done 
by the collimating force F^ - is work done against the 
magnetic pressure gradient as material is advected inwards 
near the equatorial plane. In other words it acts to increase 
the value of the Bernoulli b parameter, potentially changing 
the sign of b from negative to positive 6 . This work is then 
available to power an outflow (c.f. W01, eq. 8). 

Further study and simulation should help eluci- 
date whether turbulent magnetic pressure is indeed anti- 
collimating or collimating. It is possible the answer depends 
upon the location in the flow being examined. 

3 COMPARISON WITH COLLIMATION 
BASED ON DISPERSION RELATION 

3.1 Introduction 

As we shall see below, the estimates for the strength of the 
collimation as discussed above and in W01 and W03 are 
consistent with the estimates for the strength of the MRI 
argued by Akiyama et al. (2003). 

It is well known that fluid dynamical instabilities that 
lead to turbulence may be self-limiting, in the sense that the 
turbulence creates transport that, in many cases, damps out 
the driving terms in the instability. Thus, for example, an 
adverse entropy gradient in a stellar interior drives convec- 
tion, which in turn reduces the adverse entropy gradient. 
Similarly, it was suggested by Akiyama et al. (2003) (see 
also Akiyama & Wheeler 2003) that, under conditions in a 
core-collapse supernova where the MRI is locally unstable, 
a small seed field should grow via the MRI until the field 
reaches the critical field strength at which the MRI is no 
longer unstable. Since the instability grows exponentially, it 
was argued that it dominates linear-in-time processes such 
as passive field wrapping. The authors found that given ini- 
tial angular velocity of the supernova progenitor core of the 
order of a few radians per second, the combination of col- 
lapse and conservation of angular momentum on shells com- 
bine to create strong MRI-unstable differential rotation in 
core collapse. Further they point out that, although the ul- 
timate field orientation is not clear, a toroidal field config- 
uration seems likely, considering the large differential shear 
in core collapse. 



© 0000 RAS, MNRAS 000, 000-000 



Three routes to collimation by the MRI 9 



3.2 Estimates for the field saturation 

Akiyama et al. (2003) provide a number of different esti- 
mates for the saturation field of the MRI. One estimate 
may be obtained by setting the length scale I ~ rfi?/d(lnSl) 
equal to the characteristic mode scale. This results in a 
field strength _B 2 at ~ 4:irpR 2 Sl 2 or, equivalently, (da)^ = Vg. 
Akiyama et al. (2003) point out that this estimate of the 
field strength is substantiated by the simulations of HGB96, 
namely 27r_B s i m = B S at- In point of fact, the results of HGB96 
refer to the mean square of the field components, as the 
mean field in these simulations is identically zero, as I dis- 
cuss at length above. From the point of view of calculating 
the stresses and forces due to the field, the distinction is 
largely irrelevant. This distinction does make a large dif- 
ference, however, when it is argued that the MRI-induced 
stresses of Akiyama et al. (2003) depend upon processes such 
as helicity conservation that, through inverse cascades, cre- 
ate a significant mean field. While such processes may in- 
deed occur in core-collapse supernovae, I would argue that 
Akiyama et al. (2003) may be re-read with an understand- 
ing that the mean field may be replaced throughout their 
analysis with a tangled field with a mean stress, and this 
difference in interpretation does not change the analysis at 
all. 

Another estimate is related to this first estimate by set- 
ting the mode length equal to the local radius, resulting in 
a field strength 

B 2 ad ~ 4-KpR 2 fl{RSl' /Q.) 2 (22) 

Other estimates are found using the maximum unsta- 
ble growing mode (with a vertical seed field), or by setting 
the wavelength of the maximum growing mode equal to the 
shear length. 

Again, one of the issues that is naturally raised in their 
analysis as well as in my analysis is the question of local 
time scales and length scales. For example, the difference 
between the estimates B 2 ad and B 2 at is a factor of the ratio 
of the shear to Coriolis parameters. 

Finally, they find the ordering 

c s > rQ. ~ va- 

In particular, they find that the plasma f3 p is large, specif- 
ically the magnetic pressure is on the order of a few % of 
the gas pressure. They none the less argue that the resultant 
MHD luminosity may be significant enough to propel and 
collimate jets in core-collapse supernovae. 

Note that the similarity in the conclusions regarding the 
magnitude of the hoop-stresses based on turbulence model- 
ing and on the dispersion relation shows how the the satura- 
tion level of the MRI and the value of the effective viscosity 
is related to the MRI dispersion relation. 

Akiyama et al. (2003) do not address the magnetic pres- 
sure in their work. 



4 COMPARISON WITH SIMULATION 

Let me now briefly address recent simulations in the light of 
the preceding discussion. 

Only lately has it become possible to perform global 
DNS magnetohydrodynamic accretion flow simulations in 
three dimensions that are capable of showing the action of 



the MRI. These global simulations (Hawley, Balbus & Stone 
2001) show, as the authors point out, that the behavior of 
the flow is not adequately described by a hydrodynamic ot- 
model. This should not be surprising because, by inspection 
as I discussed above and previously, neither are local simu- 
lations adequately described by a hydrodynamic Q-model. 

Hawley, Balbus & Stone (2001) simulated the non- 
radiative magnetohydrodynamics of an initially constant 
angular momentum torus in three dimensions. They gave 
meridional contour plots of the magnetic pressure and the 
off-diagonal (R — 9) Maxwell stress, but as these plots are 
unsealed and as they do not show the other stress compo- 
nents, it is impossible to see if the MRI is creating hoop 
stresses through the tension term in the magnetic stress. 
The magnetic pressure appears to be anti-collimating (as 
I argued in W01) in the equatorial plane, and collimating 
above the plane. They found that after several orbits the 
MRI creates magnetic field with plasma /3 P on the order 
of roughly 5-10 close to the equatorial plane, whereas sev- 
eral scale heights above the equatorial plane j3 p is less than 
1. Thus, even on the equatorial plane the MRI saturates 
to a level much higher {i.e. f3 p is smaller) than in the zero 
mean-field shearing-sheet simulations of HGB96, which yield 
a plasma f3 p of roughly 40, and the MRI saturates rather to 
a level comparable to the nonzero mean-field simulations of 
HGB95, which yield a plasma f3 p of the order of 7 for the 
azimuthal field runs and of the order of 2.4 for the verti- 
cal seed field (as determined from a naive average of all the 
runs on their tables 1 and 3, and excluding run Z3). Note 
that in the case of HGB95 these values for the j3 p include 
contributions due to the mean field and the turbulent field. 

The further exploration of a very similar set of simu- 
lations in Hawley & Balbus (2002) allows us to place these 
simulations more firmly in the context of our remarks in the 
previous sections. This second set of simulations has been 
applied to the very hot (~ 10 13 K) accretion flow of Sgr A*; 
nevertheless the authors use a very simple single-component 
equation of state which allows one to estimate the ratio of 
sound speed to orbital speed from their fig. 8, namely 

c 2 ~ (1/10)« 2 . 

Using f3 p ~ 5 for the saturation of the MRI on the equatorial 
plane, then, I have 

t&~ (1/250H 2 , 

and 

2 2 2 

v A < c s < Vg . 

On the other hand, the situation changes above the equato- 
rial plane where (3 p l is much larger, with a correspondingly 
larger Alfven speed relative to sound speed. Here, v\ is likely 
a much larger fraction of Vg or vfj, but it is difficult to tell 
by how much. 

Collimation by the MRI appears to be seen in the 
2.5D simulations of Kudoh et al. (2002), Proga & Begel- 
man (2003), and Proga et al. (2003). Note that all three 
components of the velocity and magnetic field are included 
in these simulations. However, the constraint of axisymme- 
try in these simulations necessarily restricts the action of 
the MRI; in particular, this constraint eliminates the non- 
axisymmetric modes which act upon the azimuthal field, 
and which I called upon to complete the feedback loop in 
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the creation of turbulent azimuthal hoop-stresses in equi- 
librium that were the basis of my analyses. There therefore 
remains the question of the extent to which the collimation 
in Kudoh et al. (2002), Proga & Begelman (2003) and Proga 
et al. (2003) is a transient effect. This may be a moot point 
in the case of collapsars, but it is essential to the case of 
steady accretion. 

Proga & Begelman (2003) present axisymmetric sim- 
ulations of the accretion of low angular momentum gas 
onto a black hole. Despite the initial low angular momen- 
tum, at late stages of the accretion the Alfven speed, sound 
speed, azimuthal fluid velocity and Keplerian velocity are all 
roughly comparable (see their figures 9 and 10). In contrast, 
in the collapsar simulations of Proga et al. (2003), these 
velocities sometimes differ by more than an order of mag- 
nitude, and again v\ < (? s < Vg. In both cases, the radial 
pressure gradient in the equatorial plane is anti-collimating. 
Also in both cases, the field is clearly dominated by the 
toroidal component. 

Proga et al. (2003) state that it is the gradient of 
this toroidal field in particular that drives the outflow. Let 
me point out that in the meridional plane, the gradient of 
the toroidal field appears only through the magnetic pres- 
sure term. By rough examination of the printed results of 
Proga et al. (2003), it appears the toroidal field exerts a force 
through the curvature term that is comparable to the force 
due to its gradient. That is, I argue that the term B^/r in 
the divergence of the stress may be of approximately equal 
significance to the results of Proga et al. (2003) as the gra- 
dient term d r {Bi) that the authors address. (As an aside, 
note that also of comparable importance, depending upon 
location in the flow, is, by inspection, the ram pressure from 
infall, so that it appears that the jet is in part not only 
magnetically confined but inertially confined as well.) 

A final intriguing note regarding the simulations of 
Proga et al. (2003) is that the toroidal field component 
changes polarity with time and location. This confirms that 
the MRI can collimate jets through a disorganized rather 
than an organized (mean-field) toroidal field, as I suggested 
in W01 and W03. Note that the possibility of jet collimation 
by a tangled field has also been investigated in the context of 
the far-downstream jet by Li (2002) who showed that a tan- 
gled field avoids the long-mode kink instabilities that plague 
jets collimated by an ordered field. 

Finally, let me address the recent core-collapse simula- 
tions of Thompson et al. (2005) who insert MRI-induced 
viscous dissipation into ID simulations. Thompson et al. 
(2005) point out that Akiyama et al. (2003) overestimate 
the viscous time scale, and that viscous process that smooth 
out gradients in SI must be taken into consideration on the 
time scale case of core collapse. As a result, viscous heat- 
ing may significantly alter the core collapse dynamics. The 
proof-of-concept work of Akiyama et al. (2003), however, 
demonstrates that one must take into account the stream- 
wise hoop stress as well. My primary goal here and elsewhere 
has been to focus upon the hoop stress, as this had not 
been addressed in the literature. However, a full account- 
ing must at the very least address all three components, 
as I have attempted here and elsewhere, and not just one 
or another component. (The physical significance of the re- 
maining three stress components, which may come into play 
when one moves above or below the equatorial plane of sym- 



metry, is not explored here.) For example, the viscous stress 
(and not just the normal stresses) is important to the model 
of W01, as this stress transports energy and angular mo- 
mentum, which the other stresses do not. My analysis of the 
results of Thompson et al. (2005) shows that at, say, the 
fiducial radius of 20km, the inwards MRI-induced force in 
the equatorial plane, which the authors do not include in 
their code, is no more than half an order of magnitude less 
than the outwards centrifugal force that the authors do in- 
clude in their code, if one takes the modest value We = 3. 
Admittedly, both sets of simulations are severely limited by 
the inherent assumption of spherical symmetry of a ID code. 
These are valuable proof-of-principle studies; future progress 
demands multidimensional MHD codes. 



5 CONCLUSION 

In this paper I have discussed my previous work (W01, W03, 
Williams 2004a) regarding the hypothetical creation and col- 
limation of jets by MHD turbulence and by the MRI in par- 
ticular. I have added on to this work with a much more 
extensive and in-depth discussion of the dynamics of turbu- 
lent stress and the modeling of turbulent stresses in accre- 
tion flows. Also, I have placed this work within the context 
of other recent work, in particular the semianalyitic work of 
Akiyama et al. (2003) based on the dispersion relation for 
the MRI, as well as the DNS simulations of Hawley et al. 
(2001), Hawley & Balbus (2002), Proga & Begelman (2003), 
and Proga et al. (2003). 

I have focused in this paper on expounding on how sim- 
ple MHD turbulence modeling considerations in the context 
of steady shear, as discussed in Williams (2004b), lead to 
the conclusion that MHD turbulence and MRI-driven turbu- 
lence in particular can create collimating forces. As discussed 
previously (e.g. W03) I also show how collimation can be 
seen directly in the results of local shearing-sheet simulations 
of the MRI, which significantly predate global MRI simula- 
tions. Turbulent collimation by quasi-viscoelastic stresses, 
as discussed, has a natural analog in the Weissenberg ef- 
fect, whereby viscoelastic fluids climb a spinning rod. Thus, 
the turbulent MHD collimation of astrophysical jets as dis- 
cussed has a natural analog to the prehaps more prosaic 
physics of the dynamics of ordinary viscoelastic fluids such 
as egg whites which climb spinning rods. 

Indeed, this research was originally inspired by my 
serendipitous discovery in the literature of a laboratory vis- 
coelastic phenomenon well-known to rheologists (Giesekus 
(1963); Thomas & Walters (1964); National Committee for 
Fluid Mechanics Films (1972)), namely the viscoelastic flows 
discussed in W01, which seemed to bear a remarkable simi- 
larity to jet phenomena. This lead naturally to the question 
of whether turbulent conductive fluids supporting a tangled 
magnetic field can in some ways behave like viscoelastic flu- 
ids, and thus create a collimating force by tangled fields in 
analogy to the collimating force due to normal stress differ- 
ence in polymeric laboratory viscoelastic fluids. Fortuitously 
enough for this line of reasoning, simulations of the MRI in 
particular demonstrate that MRI-induced turbulent stresses 
do in fact look in many ways like the stresses in viscoelastic, 
as opposed to viscous, fluids. 

Of course, turbulent magnetized fluids are not really vis- 
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coelastic fluids, and there are many fundamental questions 
that remain regarding the behavior of MHD turbulence in 
accretion flows. Here I discuss four in particular: 

The primary open question is not whether the MRI can 
create confining magnetic hoop-stresses, but rather whether 
these stresses are strong enough to be a dominant confining 
mechanism for jet creation. Indeed, the ordering of speeds 
given in eq. (20) may be difficult to acheive in steady-state 
accretion flow scenarios. The weaker condition that the tur- 
bulent magnetic fields be dynamically significant, but per- 
haps not strong enough by themselves alone to create an 
anisotropic outflow, can be restated that the ratio Vg/v^ cw 
not be much larger than unity. Conventional wisdom would 
suggest that this condition is difficult to meet in accretion 
flows in AGN and protostellar systems. On the other hand, 
our understanding of the physics of these systems is far from 
complete. Below, I offer as a hypothesis (see hypothesis (1)) 
the notion that in these accretion systems, when jets are pro- 
duced, «t ow is a substantial fraction of, if not larger than, 
Vg, at the jet base. 

Regardless of the validity of this hypothesis, it does ap- 
pear that, as I argued in W01 and W03, the MHD turbu- 
lent hoop-stresses in accretion are at least of equal magni- 
tude to the turbulent viscous stresses. Thus, global simu- 
lations of the inner regions of accretion and jet launching 
that treat turbulence with a purely viscous a model are 
of questionable physical relevance, as pointed out by Haw- 
ley et al. (2001). For example, this result casts doubt on 
the results of, e.g., Abramowicz et al. (2002), because these 
authors do not include the turbulent hoop-stresses in their 
radial force balance equation, and as noted above, the resul- 
tant force changes the value of the Bernoulli b, which must 
be taken into consideration. Furthermore, there is a sim- 
ple and computationally cost-effective way to include these 
hoop-stresses in 2D, by the use of a viscoelastic turbulence 
model rather than a viscous model. 

A second question, from the point of view of turbu- 
lence modeling, is the relative magnitude of of the MRI 
versus other instabilities (such as convection and, possibly, 
inertial instabilities) in creating turbulent hoop-stresses. As 
has been argued previously, the MRI, as a magnetohydrody- 
namic instability rather than a purely hydrodynamic insta- 
bility, injects energy directly into the Maxwell stress tensor 
rather than indirectly, such as the case in convection, where 
energy may be expected to be transferred to the turbulent 
Maxwell stress by coupling to the Reynolds stress. Thus it is 
quite reasonable to expect that turbulent Maxwell stresses 
are relatively more significant than Reynolds stresses in 
purely MRI-driven turbulence, and this appears to be the 
case. In fact in W01 I assumed that the turbulent Maxwell 
stresses dominate the Reynolds stresses. This is not a neces- 
sary condition for the creation of collimating hoop stresses, 
but it helps, because otherwise one must contrive that the 
relaxation time for the Maxwell stresses is much longer than 
the relaxation time for the Reynolds stresses, since the ac- 
tion of the shear operator on the former creates positive 
(tension) streamwise stress whereas the action of the shear 
operator on the latter creates negative streamwise stress. In 
addition, in the case of the MRI the dimensional arguments 
presented above and in W01 and W03 argue that the ef- 



fective relaxation time is roughly equal to Q. or A. These 
arguments no longer hold when shear and Coriolis forces no 
longer drive the instability. 

A third question is the relative importance of tangled 
versus organized fields in this collimation. It remains to be 
seen the extent to which mean-field dynamo processes may 
be significant in MRI-generated collimating forces. Rather, 
here and previously I have emphasized tangled fields. I ar- 
gue collimation by MRI-induced tangled fields to be a more 
robust result than collimation by MRI-induced organized 
fields. First, I have argued above that simulations may be 
interpreted as showing that the creation of small-scale fields 
is a more reliable feature of the MRI than creation of large- 
scale mean fields. Second, if the MRI fails to create a small- 
scale tangled field, then it also fails to provide a mechanism 
for turbulent transport of angular momentum. This would 
quite clearly be a problem since it is currently accepted that 
turbulence driven by the MRI is most likely largely if not 
solely responsible for angular momentum transport in ion- 
ized disks. A third point, particularly apropos in the case 
of transient accretion events such as collapsars, is that the 
creation of a tangled field may proceed faster than the cre- 
ation of a large-scale mean field, to the extent that creation 
of a mean field depends upon inverse cascades acting over 
several orders of magnitude in length-scale. 

Thus, Lyutikov et al. (2003), in the context of gamma- 
ray burst (GRB) fireballs, argue that the prompt creation of 
a large-scale field is problematic. However, citing Gruzinov 
& Waxman (1999) who in turn cite Loeb & Perna (1998), 
they argue that the linear polarization of GRB0212206 seen 
by the RHESSI satellite as described by Coburn & Boggs 
(2003) implies a large-scale field [as indeed Coburn & Boggs 
(2003) argue as well] and that, since it is difficult to create 
such a field quickly, the field must exist prior to the catas- 
trophe producing the GRB: 'Such fields cannot be generated 
in a hydrodynamically-dominated outflow, which is causally 
disconnected on large scales. Thus, the large scale magnetic 
fields should be present [at the beginning].', Lyutikov et al. 
(2003). The authors claim that '[A large polarization] im- 
plies that magnetic field coherence scale is larger than the 
size of the visible emitting region. . . '. The truth of this as- 
sertion depends upon what is meant by 'coherence scale'. In 
particular, it is wrong to assert that the polarization implies 
a large-scale field: tangledness should in no way be confused 
with isotropy. 

As an analogy, the orientation of mesogens in the 
nematic phase of liquid-crystal displays (LCDs) is coher- 
ent on a macroscopic scale. This large-scale anisotropy (or 
aeolotropy) makes the speed of light depend upon polar- 
ization, even though the size of an individual mesogen is 
typically smaller than an optical wavelength (Rosenblatt, 
private communication). Likewise, a large-scale mean dyad 
field BtBj does not imply a large-scale mean vector field 
Bi. Processes that do not discriminate between B and — B 
care about the former, not the latter; this includes the lin- 
ear polarized emission of cyclotron and synchrotron radia- 
tion. In the case of cyclotron and synchrotron radiation, only 
the presence of circular polarization discriminates between 
a mean field and a mean dyad (Matsumiya & Ioka (2003)). 
Thus, large linear synchrotron polarization is possible due to 
a tangled field, so long as Mij is strongly anisotropic (Sagiv 
et al. (2004)). 
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Again regarding tangled versus ordered fields, Thomp- 
son et al. (2005) state, regarding the work of Akiyama et al. 
(2003), that '[The argument presented by Akiyama et al. 
(2003)] assumes that the magnetic field generated by the 
MR1 can form the organized large-scale fields required for 
collimation and jet formation.' Strictly speaking this may 
be true, since Akiyama et al. (2003) refer to mean fields 
rather than tangled fields. But, inasmuch as Akiyama et al. 
(2003) may be reinterpreted as describing the magnitude of 
the stress due to a tangled field rather then an ordered field 
as described above, the estimates for the forces due to the 
MR1 remain the same. Thus, while the authors of Akiyama 
et al. (2003) were perhaps unaware of it, the processes they 
described does not actually depend upon the creation of or- 
ganized large-scale fields, and in this sense Thompson et al. 
(2005) are incorrect. 

Fourth and finally, the role of magnetic turbulent pres- 
sure and other stresses such as the neglected inertial stresses 
requires further analysis. I have argued here that the turbu- 
lent magnetic pressure will be anti-collimating in jet scenar- 
ios, and that it provides a reservoir of energy to accelerate 
the jet, but this conclusion is more tenuous than the more 
robust conclusions regarding the collimating normal stress 
differences caused by MHD turbulence. In addition, is con- 
ceivable that neither magnetic stresses nor inertial stresses 
are large enough by themselves to launch and collimate jets, 
but that the two acting together may create jets. Discus- 
sion of intertial stresses in particular is left to future work; 
notwithstanding their potential dynamical significance, I hy- 
pothesise below that the MHD turbulent stresses must be 
dynamically significant in a broad range of jet-producing 
accretion flows. 

These problems discussed above, as well as many oth- 
ers, remain to be answered completely. None the less, I have 
shown here that my previous estimates of the turbulent colli- 
mating force are consistent with other estimates, and in par- 
ticular the three lines of evidence and reasoning I presented 
above all reduce to the same condition, which may be stated 
in two different ways: The energy in the toroidal field is a 
sizeable fraction the orbital kinetic energy, or equivalently, 
the local toroidal magnetic wave speed is a sizeable fraction 
of the orbital speed. The difference between the theoreti- 
cal approach based on the MRI dispersion relation and the 
approach based on turbulence modeling is simply whether 
the wave in question is a purely Alfven wave or an Alfven- 
like turbulent transverse elastic wave. It is not yet clear how 
sizeable a field is necessary, and whether condition (20) need 
indeed be met anywhere in the flow for a turbulent toroidal 
magnetic field to be dynamically significant. 

Notwithstanding these remaining questions, then, I of- 
fer as hypotheses the notions that: 

(i) MHD turbulence generated by the MRI (as well as 
perhaps additional instabilities) creates dynamically signifi- 
cant toroidal hoop stresses (more properly, a toroidal normal 
stress difference) that contribute significantly to jet creation 
and collimation; 

(ii) this process does not depend upon mean-field dy- 
namo processes, and in particular much of the field energy 
may be in a tangled field rather than a mean field; 

(iii) inasmuch as MHD turbulence (driven by the MRI 



in particular) is thought to be the predominant source of 
turbulent viscous angular momentum transport in hot ion- 
ized disks and accretion flows, the corresponding turbulent 
elastic collimation process described here may be a nearly 
universal aspect of jet formation, applicable not just to the 
transient dynamics of collapsars but to jets in a broad range 
of quasi-steady- state systems such as AGN and quasars, 
microquasars, and protostellar systems; this hypothesis is 
offered notwithstanding the predominant understanding of 
thick accretion flows that would argue that such azimuthal 
turbulent magnetic fields should be relatively weak; 

(iv) that this jet creation and collimation, as described, 
is not offered as a process that occurs in addition to tra- 
ditional magnetocentrifugal acceleration mechanisms, but 
rather as an alternative to them. 

Blandford & Payne (1982), Shu et al. (1994) and oth- 
ers have argued that jets are created by material flung out 
following large-scale open magnetic field lines. To be quite 
explicit, as a hypothesis, I reject the notion that such mag- 
netocentrifugal mechanisms are responsible for tightly colli- 
mated jets. In such scenarios, power can not be transferred 
to the jet without transferring angular momentum as well. 
Note that recent observations of certain jets suggest that 
the fields may have helical structure (Gabuzda, Murray & 
Cronin (2004)) and the jets may be rotating (Coffey et al. 
(2004)). However, this in itself does not demonstrate mag- 
netocentrifugal acceleration. Indeed, it would be surprising 
if Nature conspired to remove all of the residual angular mo- 
mentum from material before accelerating it in a jet. The 
presence of a helical magnetic field, rather than a purely az- 
imuthal or vertical field, can be explained as the result of 
advection, and does not necessarily say what is the dynami- 
cal role of that field, and in particular it does not necessarily 
imply that material is in some sense being flung out along 
field lines. Rather it is suggested here that the predominant 
role of the magnetic field in the central engine of a broad 
array of jet sources is the creation of turbulent hoop-stresses 
as described, which create and collimate an outflow deep in 
the core of the central engine. In analogy to viscoelastic lab- 
oratory phenomena, and in contrast to magnetocentrifugal 
mechanisms, in such a scenario material can take part in a 
collimated outflow not because angular momentum is trans- 
ferred to it, but because angular momentum is removed from 
it. 
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